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ABSTRACT 

We measure shifts of the acoustic scale due to nonlinear growth and redshift distortions to a 
high precision using a very large volume of high-force-resolution simulations. We compare results 
from various sets of simulations that differ in their force, volume, and mass resolution. We find a 
consistency within 1.5 — cr for shift values from different simulations and derive shift oi{z) — 1 = 
(0.300 ± 0.015)%[Z?(z)/£'(0)]^ using our fiducial set. We find a strong correlation with a non-unity 
slope between shifts in real space and in redshift space and a weak correlation between the initial 
redshift and low redshift. Density-field reconstruction not only removes the mean shifts and reduces 
errors on the mean, but also tightens the correlations. After reconstruction, we recover a slope of 
near unity for the correlation between the real and redshift space and restore a strong correlation 
between the initial and the low redshifts. We derive propagators and mode-coupling terms from our 
A^body simulations and compare with the Zel'dovich approximation and the shifts measured from 
the fitting, respectively. We interpret the propagator and the mode-coupling term of a nonlinear 
density field in the context of an average and a dispersion of its complex Fourier coefficients relative 
to those of the linear density field; from these two terms, we derive a signal-to-noise ratio of the 
acoustic peak measurement. We attempt to improve our reconstruction method by implementing 
2LPT and iterative operations, but we obtain little improvement. The Fisher matrix estimates of 
uncertaint y in the acous tic scale is tested using 5000/i~'^ Gpc'^ of cosmological PM simulations from 
[Takahashi et al.l (|2009a[ ). At an expected sample varianc e level of 1%, the agreement between the 
Fisher matrix estimates based on iSeo fc EisensteinI (|2Q07[ ) and the A-body results is better than 10 
%. 

Subject headings: distance scale- large-scale structure of universe - methods: numerical 



1. INTRODUCTION 

In recent years, attention to baryon acoustic oscil- 
lations (BAO) as a dark energy probe has increased 
unprecedentcdly due to its robust nature against sys- 
tematics, and they are now an essential component of 
most of the major future dark energy surveys under 
consideration. BAO originate from the sound waves 
that propagated through the hot plasma of photons and 
baryons in the very early Universe. At the epoch of 
recombination, photons and baryons decouple, and as a 
result, the sound waves freeze out, leaving a distinctive 
oscillatory feature in the large-scal e structure of the 
cosmic microwave backg round fe.g.. iMiller et al.l 119991: 
de Bernardis et all 120001 : iHananv et al.l 200^ Lee et al.l 



20011: IHalverson et all I2002t INetterfield et al.l I2002l:_ 
Pearson c t all 120031: iBenoit et all 120031 : iBennett et al.l 
2003: Hinshaw et al.l 12007. JOOj) and the matte r den- 
sity fields in Fourier space fe.g.. [Peebles k, Yvi 



Sunvaev fc ZeldovichI [l970l: iBond fc Efstathiou 



1970; 



1984; 



HoltzmanI 119891 : IHu fc Sugivamal Il996t IHu fc Whit^ 
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199e : lEisenstein fc Hul[T99llMeiksin. White, fc Peacockl 
1999[ ): In configuration space, the BAO appears as a 
single spherical peak at its characteristic scale. The 
characteristic physical scale of this oscillatory feature, 
BAO, is the distance that the sound waves have traveled 
before the epoch of recombination, which is known 
as "sound horizon scale" . This sound horizon scale 
is and will be measured precisely from current and 
future CMB data. With the knowledge of the physical 
scale, BAO can be used as a standard ruler to measure 
angular diameter distance and Hubble parameter at 
various redshifts and therefore provi de critical informa- 
tion to id e ntify d ark energy fe.g.. IHu fc White! 119961: 



EisensteinI [20031: [Blake fcGlazebrookl l2003t iLindeH 
20031: IHu fc HaimanI 120031 ISeo fc EisensteinI I2003D . 
Recently, BAO have been detected from large-scale 
structure of galaxy distributions and have been used 



(lEisenstein et al.l [2005: Cole et al. 


120051 IHiitsil 120061: 


Tegmark et al.l 20061 IPercival et al. 


2007allbl: Blake et al. 


2007 


; Padmanabhan et al.l 2007; Okumura et al. 


2008 


; [Estrada et all 20091: iGaztaiiaga et al. 200911?: 


Sanchez et all 120091: IPercival et al.l 120101: iKazin et al. 


20101). 



Due to nonlinear structure growth at late times, the 
oscillatory feature of the BAO is increasingly damped, 
proceeding from small sca les to larger scales, with 
decreasing redshift fe . g., iMeiksin. White, fc Peacockl 
19991: ISeo fc EisensteinI i2005t iJeong fc Komat su 200e 



Eiscnstcin et al. 120071 : ICrocce fc Scoccimarrg |200 



Matsubaral I2008f) . Redshift distortions enhance 
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nonlinear damping along th e l ine of si ght direction(e.g., 
Meiksin, White, fc Peacockl 119991: iSeo fc EisensteinI 



200a lEisenstein et al.ll2'007t iMatsubaral I2008D . Despite 



the resulting loss of the BAO signal with decreasing 
redshift, BAO are believed to be a robust standard 
ruler. The sound horizon scale is well determined from 
the CMB and the scale corresponds to ~ lOO/i^^ Mpc 
in present time, implying that the feature is still mostly 
on linear scales where evolutionary and observational 
effects are much simpler to predict than in the nonlinear 
regime. Indeed, such degradation in the contrast of 
the BAO due to nonlinear structure growth, redshift 
distortions, and possibly galaxy bias has been studied 
since the late 90's and is relati v ely well understood 
(iMeiksin. White, fc Peacock [l999l: iSpringel et all [20051: 



Angulo et al 



2005t I Eisenstein et al 



ISeo fc EisensteinI boOsF 



Crocce fc Scoccimarrc 



iJeong fc Komatsj] 



Eisenstein et al.l 




Huff et al.ll2007HSmith"c t al. 2057l: iMatarrese fc Pietronil 



20071: iNisMmichi et ai.i ,20071: iSmith et al I 



Angulo et all 120081 (C rocce fc Scoccimarrol 



2008 



2008 



Matsubar al 120081: iT akahas hi et al.l 120081: 'Sanche z et alj 
2008; Jco ng fc Kom atsu 20091: iTaruva et a l. 2009). 

Nonlincarity also induces a shift of the BAO scale in 
the low-redshift matter distribution relative to the sound 
horizon scale measured from the CMB (e.g.. ISmith et al.l 
l2008t ICrocce fc Scoccimarrol [200l [Sanchez et al.ll2008D . 
As the demand for acoustic peak accuracy moves from 
the current level of a few % to the sub-percent level of 
future surveys, we are now required to understand the 
systematics on the BAO to much better precision. While 
it is evident that the shift of the BAO scale depends on 
the choice of the estimator, most recent studies seem 
to lay more weight on residual shifts using optimal es- 
timators being at a sub-percent level at z ^ when 
accounting for nonlinear structure grow t h and redshift 
disto r tions (ICrocce fc Scoccimarrol 120081 : [Sanchez et al.l 
[200a ISeo et al.l I2008D. or eve n with halo/galaxy bias 
([Padmanabhan fc Whitil2009D . Certainly more studies 
are necessary to confirm the results and ultimately reach 
a general consensus. 

In the previous study (|Seo et al.l l2008f ) (hereafter 
SSEW08), we investigated effects of nonlinear evolution 
on BAO using a large volume (i.e., 320h~^ Gpc^) of PM 
simulations. We found that the shift on the acoustic scale 
indeed increases with decreasing redshift and is less than 
a percent even at z = 0.3 in redshift space. 

In this work, we extend the study by using high-force- 
resolution simulations that are generated by a new N- 
body code ABACUS by Metchnik fc Pinto (in prepara- 
tion). With the new simulations, we update the evo- 
lution of the shifts on the acoustic scale due to nonlin- 
ear structure growth and redshift distortions for various 
force, volume, and mass resolutions. We calculate prop- 
agators, i.e., the correlations between the linear density 
fields and the nonlinear density fields at low redshift, 
which directly manifest the nonlinear damping of the 
BAO. We also derive the mode-coupling contribution to 
the power spectrum with an attempt to qualitatively re- 
late this to the measured v alues of the shifts. Recently, 
[Padmanabhan et al.l (|2009( ) showed that the density field 
after reconstruction is not the linear density field at sec- 
ond order. We relate the propagator and the mode- 



coupling term to an average and a dispersion of non- 
linear density fields in the complex Fourier plane before 
and after reconstruction, relative to the linear density 
fields, and derive a signal-to-noise ratio of the standard 
ruler test from these two terms. The effect of galaxy bias 
will be prese nted in compan ion papers (Mchta et al. in 
preparation; IXu et al.]|20100 . 

It has been demonstrated that the original density-field 
reconstruction scheme based on the Ze l' dovic h approx- 
imation, presented by lEis enstein et al.l ([20071 1 , is quite 
efficient for removing nonlinear degradation o n BAO de- 
spite its simplicity ( Seo fc EisensteinI \200% iHuff et al.l 



[20^ [Seoetall 120081 IPadmanabhan et al.l I2009D . Efln 
ciency of reconstruction in terms of an increase in the 
signal-to-noisc depends on the redshift and the shot noise 
of the density fields. Meanwhile, it has been shown that 
the reconstruction removes almost all of the nonlinear 
shifts of the acoustic scale even when it seemingly is not 
efficient in terms of the signal-to-noise ratio (SSEW08). 
This success, conversely, can be interpreted that a fur- 
ther improvement on the reconstruction scheme will be 
only a second order effect, at least for BAO. Neverthe- 
less, we discuss and test possible improvements on our 
fiducial scheme. While there are other sophisticated re- 
construction methods in the literature aimed for the re- 
covery of velocity field s and the initial density fields (e.g., 
iMohavaee et al.H2006( ). in this paper we limit ourselves to 
mild modifications to our fiducial method, mainly due to 
its proven success for BAO. We first test an implcmcn- 
tation of 2LPT in stead of the Zel'dovich approximation 
([Zel'dovich! [19700 and second, test iterative operations 
of the fiducial method in order to improve the signal-to- 
noise ratio. 

The importance of an accurate prediction of the signal- 
to-noisc ratio for future BAO surveys is evident. The cal- 
ibration of the Fisher matrix-based estimations against 
the A^-body results have been tried repeatedly, and the re- 
sulting discrepancy is at most 20% (e.g., SSEW08). Fur- 
ther calibration is often limited by the volume of the sim- 
ulations: in order to minimize the dispersion of the dis- 
persion among different measurements, we need a large 
number of random subsamplcs while each subsample has 
an enough cosmic volume to measure the BAO scale. 
In this pa per, we calibrate t he Fis her matrix estimation 
based on ISeo fc EisensteinI ([20071) to a level of 1% by 
utilizing t he enormou s cosmi c volume of 5000h~^ Gpc"^ 
from Takahashi et all ([2009a[) . 

This paper is organized as following. In § [2l we de- 
scribe our new cosmological iV-body simulations and the 
methods of analysis to measure the acoustic scales 
from the simulations. In § [3l we present the resulting 
shifts and errors on the measurements of the acoustic 
scale when accounting for the nonlinear growth and red- 
shift distortions before and after reconstruction. In § [4| 
and § [SI we derive propagators and mode coupling terms 
from the simulations and qualitatively compare these 
with the errors and the mean values of the shifts from 
the simulations. In § [SI we relate the propagator and 
the mode-coupling term to the signal-to-noise ratio of 
the standard ruler test. In § [7l we implement 2LPT and 
iteration into our reconstruction scheme and discuss the 
eff ect. In [HI we use the SOOO/i""^ Gpc"^ of simulations 
bv lTakahashi et all ([2009a[ ) to test the Fisher matrix cal- 
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culations in iSeo fc EisensteiiJ (|2007f ) given the minimal 
sample variance. Finally, in §[9l we summarize the major 
results presented in this paper. 

2. METHODS 
2.1. Simulations 

We use three sets of iV-body simulations to produce 
the results presented in this paper. The first two are 
produced using a new high-force resolution iV-body code 
ABACUS by Metchnik & Pinto (in preparation) and 
are used for measuring the shifts of the acoustic scale 
at different redshifts. Our ABACUS simulations have 
high force resolution, sufficient to resolve dark matter 
halos, unlike the particle-mesh simulations in SSEW08. 
The ABACUS code uses a new method to solve the far- 
field gravitational force, resulting in higher force accu- 
racy than standard tree and Fourier methods, at high 
computational speed: the code speed is such that each 
1024^ particle simulation takes only 3 days on a single 
16-core machine. 

We generate these two sets using cosmological param- 
eters that are consistent w ith WMAP5+SN+BA0 re- 
sults (jKomatsu et al.ll2009[) : 0^ = 0.279, Qa = 0.721, 
h = 0.701, VLb = 0.0462, = 0.96, and CTg = 0.817. The 
initial conditions are generate d at z = 50 using the 2LPT 
based IC code bv lSnkol (l2005l) : we do not correct for the 
finite volume effect on the DC mode but assume a zero 
DC mode. The two sets differ in their volume and mass 
resolution. 

The first set employs 576'^ particles in a box of 
2/i~i Gpc and a particle mass of 3.2413 x 10^^ hr'^MQ, 
with a softening length of 0.1736/1"^ Mpc (i.e., 1/20 of 
the interparticle spacing) for gravity calculation. We 
generate a total of 63 boxes and therefore a volume of 
504/i~^ Gpc^. We compute power spectra at z = 3.0, 
2, 1.5, 1, 0.7, and 0.3 using 576^ density grids using 
the cloud-in-cell interpolation. The power spectra are 
spherically averaged within wavenumber bins of width 
Afc = 0.001ft. Mpc~^ in real space and redshift space. 
We will refer to this sample as "G576" . 

The second set employs 1024^ particles in a box of 
\hr^ Gpc and a particle mass of 7.2110 x 10^°/i"1Mq, 
with a softening length of 0.0488/i~^ Mpc for gravity cal- 
culation. We generate a total of 44 boxes and there- 
fore AAh~^ Gpc^. The power spectra are computed 
in the same way as for G576 at z = 1 (i.e., Afe = 
O.OOl/i Mpc~^), but using 1024"^ density grids this time. 
We will refer to this sample as "G1024". We use G1024 
to test whether any of the results depend on the different 
force or mass resolution or simulation volume. 

For t he thir d set, we use the A^-body results generated 
bv lTakahashi et al.l (|2009a[). This simulation is generated 
using the Gadget-2 code (jSpringel et al.|[200lt iSpringell 
[2005[i in a Particle-Mesh (PM) mode wi th cosmologi- 
cal p arameters based on WMAP3 results (jSpergel et al.l 
[2003): = 0.238, fl^ = 0.762, h = 0.732, = 0.041, 
Us = 0.958, and as = 0.76. The set employs 256"^ par- 
ticles in a volume of lh~^ Gpc, and therefore a particle 
mass of 3.94 x 10^^ h~^MQ, and a 256^ force mesh. The 
initial condition is generated at z = 20 using Zel'dovich 
approximations, and power spectra are computed at 
z = 3, 1, and using 512^ grids and are spherically 



averaged for each discrete value of radial f^. the width 
Afc = 0. While this set has worse force resolution, the 
virtue of this sample is its large volume: 5000h~^ Gpc'^! 
We will refer to this sample as "T256". We utilize this 
set, first to check the consistency with respect to G576 
and G1024, but more importantly to calibrate the Fisher 
matrix error estimates, that is, to measure to high preci- 
sion the scatter of the acoustic scale shift. We summarize 
our three A^-body sets in Table [1] 

2.2. Fitting methods 

To measure the shift of the acoustic scale, we fit mod- 
els to the spherically averaged power spectra both in real 
and redshift space. The fitting method used in this pa- 
per is identical to SSEW08. To summarize, we fit the 
observed power spectra Pobs to the following fitting for- 
mula: 

Pobs(fc) = B{k)Pm{k/a) + A{k), (1) 

where a, B{k), and A{k) are fitting parameters. Here a 
is a scale dilation parameter and represents the ratio of 
the true (or linear) acoustic scale to the measured scale. 
For example, a > 1 means that the measured BAO be- 
ing shifted toward larger k relative to the linear power 
spectrum. If the nonlinear shift of the acoustic scale 
were not corrected, a would represent the ratio of the 
mis-measured distance to the true distance. The tem- 
plate power spectrum Pm is generated by modifying the 
BAO portion of the linear power spectrum with a nonlin- 
ear parameter Sm to account for the degradation of the 
BAO due to nonlinear growth and redshift distortions as 
discussed in SSEW08. 

Praik) = [Fii„(fc) - F„w(fc)] exp [-k^^lJ2] 

+P„w(fc), (2) 

where Pun is the li near power spectrum an d Pnw is the 
nowiggle form from lEisenstein fc Hul (|1998| ). We do not 
include S free parameter. As demonstrated in 

SSEW08, our results are not sensitive for a wide range 
of !],„ (i.e., AEm = ±2h-^ Mpc for z < 3.0) due to 
the large number of free fitting parameters allowed in 
B{k) and A{k). We therefore simply choose fiducial val- 
ues of Sm based on th e Zel'dovich approximation from 
lEisensteiii et al.l (|2007| ). 

The term B{k) allows a scale-dependent nonlinear 
growth, and A{k) represents an anomalous power, i.e., 
additive terms from the nonlinear growth and shot noise. 
By including both B(k) and A{k) with a large num- 
ber of parameters, we minimize the contribution to the 
standard ruler method from the broadband shape of the 
power spectrum. We adopt three different choices of 
parametrization for B{k) and A{k). For real space, we 
use the following two choices. For the first choice, which 
we call "Poly 7" , we use a quadratic polynomial in k for 
B{k) and a 7*'' order polynomial function for A{k). For 
the second choice, referred to as "Pade", we use Pade 
approximants for B(k) in a form of 6o(l -I- Cik + c^k^ + 
c^k^) /(I + C2fc 4- c^k'^), and a quadratic polynomial in k 
for A{k). 

^ The wavenumber is discretized as fc = 27r/Lijox x (n^ -|- n| + 
n^)^^^ where Lbox = l/i^^ Gpc and {i=l-3) is an integer. 
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Fig. 1. — Power spectra divided by a smooth, nowiggle power spectrum Pnw at z=3.0, 1.0, and 0.3 in real (left) and redshift space{right), 
before (black for G576 and red for G1024) and after reconstruction (b lue for G576 and m agenta for G1024). The dashed line is the linear 
theory model. The dotted lino at z = 1 is from the halofit model from lSmith et aTl I I2003I) . We do not subtract shot noise. 



For the redshift space fits, we extend the parametriza- 
tion in equation ([T]) to include a finger-of-God (FoG) 
term. 

Pik) - [B(A:)P„ (k/a) + Aik)] x Ffog (3) 
+ ei, 

with Ffog = cxp [-{kdif^] . 

(4) 

Here B{k) and A{k) are quadratic polynomials in k, 
and di > Oh~^ Mpc and ^2 > 0. We caU this "Poly- 
Exp-Out". Wc use a fitting range of 0.02/i Mpc~^ < 
k < 0.4/iMpc~^. Our results are highly consistent 
for 0.35/i Mpc"^ < fcmax < 0.5/iMpc"i and fc,nin < 
O.OAh Mpc^^ where /cmin and /cj^ax sue the two outer 
bounds of the fitting range. 

As demonstrated in SSEW08, the fits are virtually the 
same for different choices of fitting formulae given the 
sample variance. While we try variations in Em and fit- 
ting formulae for all the power spectra we generated, we 
only quote values derived using Poly? in real space and 
using Poly-Exp-Out in redshift space. 

2.3. Resampling methods 

The true covariancc matrix of the A'^body simulations 
is unknown a priori. Therefore we use the variation be- 



tween the simulations to assess the true scatter in a. We 
fit each simulation assuming that the covariance matrix 
for P{k) is that of a Gaussian random field, i.e., assuming 
independent band powers with variances determined by 
the number of independent modes in each. Although this 
is not an optimal weighting of the data for the determina- 
tion of a, the effects from non-Gaussianity in the density 
field will still be reflected in the scatter between best-fit 
g's fr om different simulations (also see iTakahashi et al.l 
120091 for the negligible effect of non- Gaussian errors in 
multi-parameter fitting). 

We use various resampling methods to measure the 
shifts and the scatters of shifts of the acoustic scale. In 
the first method, wc generate 1000 subsamples by ran- 
domly resampling M boxes out of total N simulations 
without replacement. We perform a analysis to the 
individual subsamples and find the mean and the scat- 
ter in the best fit a. The scatter in a is rescaled by 



V M /y/N — M to give a scatter associated with the mean 
value of a for a total of N simulations, and is presented 
in Table EI We choose M = 31 for G576 {N = 63) and 
M = 22 for G1024 {N = 44). Note that M ~ N/2, and 
now the scatter of the subsamples (presented in Figure 
[21 m and El) is very similar to the scatter associated with 
the mean value of a (in Table ^ . 
We also use jackknife resampling to measure the shifts 
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and find that both methods give consistent estimates 
of the shift and the scatter of the shift for G576 and 
G1024. For T256, we measure the shifts primarily by 
using the jaekknife resamphng; we also use a third re- 
sampling method, Bootstrap resampling, to measure the 
scatters among 125 subsamples of 40 boxes and 250 sub- 
samples of 20 boxes. For T256, the means of the shifts 
from the two resampling methods are almost identical 
and the errors on the shifts agree within their expected 
uncertainties. 

2.4. Reconstruction method 

We reduce degradations in the BAO due to nonlin- 
ear growth and redshift distortions us ing the density- 
field r econstruction method described in lEisenstein et ahl 
(|2007[ ). We smooth gravity of the observed nonlinear 
density fields on small scales with a Gaussian filter, 
apply a linear perturbation theory continuity equation 
V ■ q ~ —S, and derive the linear theory motions q. We 
displace the real particles and the uniformly distributed 
reference particles by ~q, derive the two density fields, 
and find the final, reconstructed density field by sub- 
tracting the two. We use R = 10h~^ Mpc as our fidu- 
cial width of the Gaussian filter. In redshift space, the 
displacement field derived from the observed redshift- 
space density field is multiplied by (1 -I- /)/(! -I- /3) along 
the line of sight in order to approximately account for 
the linear bias and the redshift distortions. For mass, 
(1 + /)/(! + (3) = 1, meaning that the redshift-space 
density field exa ctly account s for the redshift distortions 
in linear theory (|Kaiserlll987f n. We do not correct for the 
FoG compression, but the effect of the FoG compression 
is briefiy discussed in § [S] 

3. THE EVOLUTION OF ACOUSTIC SCALES 

We show the spherically averaged real-space and 
redshift-space power spectra of the matter distribution 
for G576 and G1024 in Figure [T] As expected from its 
higher mass resolution, G1024 (at z = 1) shows slightly 
higher power on small scales than G576 in real space 
while showing more FoG effect in redshift space. 

Figure [2] and Table [2] show the resulting mean and the 
errors of the mean of the shifts from the subsamples, in 
comparison to the results from T256 and previous results 
from SSEW08. We find that the shifts measured from 
G576 are less than those measured in SSEW08 by up to 
~ 2.7 — cr. However, the cosmologies in the two sets are 
not identical, which we need to take into account. We de- 
rive a cosmological scaling prediction for shifts using the 
se cond order perturbation theory (hereafter 2PT) based 
on iPadmanabhan &: Whit^ (|2009| 1 and compare the pre- 
dictions (broken lines with the corresponding color in the 
figures) with the measurements. Indeed the shifts pre- 
dicted for the cosmology of SSEW08 are higher than for 
G576. Then the disagreement between SSEW08 and the 
prediction is about 1.8 — a. 

The measured shifts for G576 and T256 are in a good 
agreement (i.e., within 1 — cr at the shared redshift bins) 
and the 2PT prediction matches the numerical result to 
a reasonable level, especially for z > 1. The agreement 
appears worse at low redshift, which might be relevant 

Here / is the logarithmic derivatives of the linear growth factor 
and /3 = f /bias. 



to 2PT being a better approximation at high redshift. 
Meanwhile, the shift aX z — 1 from G1024 is 1.5 — cr 
smaller than G576. In redshift space, the shift values 
from G1024 are bigger than that in real space, as ex- 
pected, and are in a better agreement with other sets 
than in real space (i.e., within 1 — cr). As the sample vari- 
ance of G1024 is relatively large compared to the other 
sets, we suspect that the discrepancy is largely due to the 
statistical fluctuations rather than a resolution or simu- 
lation volume. 

The fig ure illustrates th a t our reconstruction scheme 
based on lEisenstein et al.l (|2007| ) reduces the shifts es- 
sentially to zero in most of the cases, while the shift is 
consistent with zero at a 1 — cr level for G1024 in real 
space. The errors on shifts have decreased by a factor 
of 2 (2.6) at z = 0.3, a factor of 1.5 (2) at z = 1, and 
a factor of 1.2 (1.4) at z = 3 in real space (in redshift 
space). 

Left panels of Figure [3] and 2] show the distributions 
of a — 1 for all the subsamples for G576 before recon- 
struction. An intriguing feature is the strong correla- 
tion between a — 1 in real space and a — 1 in redshift 
space in the lower left corners of the figures. The slope 
of the correlation is larger than unity and this implies 
that shifts in redshift space are greater than shifts in 
real space. The strong correlation and the non-unity 
slope is expected, as the displacements in redshift space 
are correlated with bu t larger than those in real space 
(jEisenstein et al.ll2007f ). Meanwhile, there seems no ob- 
vious correlation of a — 1 between the low redshift and 
the initial redshift. While not shown in the figure, the 
correlation between any two redshifts becomes stronger 
for a smaller redshift difference. 

With reconstruction, the right panels of the figures 
show that we recover a slope of near unity for the 
correlation between the real and redshift space a — 
1. Also note that we recover a strong correlation be- 
tween the low and the initial redshifts. The non- 
linear growth involves a second-order process that is 
imperfectly correlated with the initial density fields, 
and the former dominat es over the shifts imprinted 
by the initial condition (iCrocce fc ScoccimaTTol 120081 : 
IPadmanabhan k, Whitell2009() . As the reconstruction re- 
moves shifts due to this second-order process, i.e., due 
to nonlinear structure growth, the remaining shifts are 
dominated by the variations in the initial conditions. 

Because the various a's within individual subsamples 
are highly corrected, it is effective to study the differ- 
ences of a's, as the errors on the differences is consider- 
ably smaller than the quadrature sum of the two terms. 
This allows us to study certain trends with reduced sam- 
ple variance. We derive the differences of a, that is, 
Aa^r = Cis — Cir, whcrc as is for redshift space and 
for real space, and Aaz^zso = o-z — oizso, where az = otr 
or a^, for each subsample and calculate the mean and 
the errors on the mean. The results are listed in Ta- 
ble [2] and are also graphically represented with the gray 
diagonal lines in Figure [3] and |4l the y-intercept or x- 
intercept of the middle gray line shows Aa and the outer 
gray lines depict a 1 — cr range for Aa. The results in 
Table [5] are consistent with the trends of correlation we 
observed in the figures: errors on Aa^r are indeed less 
than a quadratic sum of and a^, implying a strong 
correlation between shifts in real and redshift space. Er- 
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Fig. 2. — The nonlinear evolution of shifts a — 1 with redshift. Top panels show the mean and the errors of the mean of a before 
reconstructi on. The broken lines with the co rresponding colors arc the expected nonlinear shifts derived using the second order perturbation 
theory, as in lPadmanabhan fc Whitd l|2009l ). In the bottom panels, solid points show the values after reconstruction, in comparison to the 
values before reconstruction (open points). Data points for G1024 arc slightly displaced in z for clarification. We note that the sample 
variance is highly correlated between shifts at different redshifts for a given set of simulations. The solid red line in the bottom panel is 
our best fit to a(z) — 1 oc [D{z)/ D{0)]^ when including the covariance between redshifts. 
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rors on Ackz^zso before reconstruction, on the other hand, 
are very similar to a quadratic sum, which is consistent 
with the httle correlation observed in the figures. Af- 
ter reconstruction, errors on all differences become much 
smaller than the quadratic sums. 

Figure [5] shows results for G1024 (black points and red 
error bars) that are very similar to G576 (blue error bar) 
except that the mean of the distribution seems biased rel- 
ative to G576 by 1.5 — tr. We believe that sample variance 
can be a reasonable cause for 1.5a discrepancy between 
G576 and G1024. Obviously, reconstruction cannot cor- 
rect shifts in this case: note that the sample variance 
on shifts after reconstruction is at the level of the initial 
condition, and reconstruction cannot do anything to the 
sample variance in the initial conditions. 

Finally, we attempt to model the growth of a with 
redshift as a function of the linear growth factor D{z). 
In order to reduce the sample variance effect from the 
initial condition, we use Aa^^^so to derive the growth 
of a. From the results of G576, when we perform a 
analysis to the mean a at z = 0.3, 1, 3.0 with the 
covariance derived from the resampling, we find that 
a{z) - 1 = (0.295 ± 0.075)%[D(z)/D(0)]i-74±o.35^ The 
resulting power index is in agreement with the expected 
D'^{z) dependence from the 2P T within 1 — a, and is in 
excelle nt agreement with what iPadmanabhan fc White! 
(|2009f ) measured using their fitting, despite the dif- 
ferent cosmology and different simulation sets used. It 
is also in good agreement with the power index mea- 
sured from SSEW08. If we fix the power index to be 
2, as expected from the 2PT, we find the best fit of 
a{z) - 1 = (0.300 ± 0.015)%[D(z)/D(0)]2 (solid line in 
Figure [11). 

Figure [5] shows a in real space as a function of D'^{z) 
by which we absorb differences due to different growth 
factors between cosmologies. The 2PT predictions ap- 
pear in straight lines in this case and the different slopes 
are due to the different nonlinear structure growth (i.e., 
effects of the second order term on the BAO) depending 
on cosmology. Now the predictions for T256 and G576 



1 — ' — ' — ' — ' — I — ' — ^ 

o G576, 504h-='Gpc3 

G1024. 44h-3Gpc=i 

X SSEW, 320 
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Fig. 6. — We replot the top left panel of Figure [2] as a func- 
tion of D'^{z), where D{z) is the linear growth factor, in order to 
absorb differences due to different growth factors between cosmolo- 
gies. The 2PT predictions (broken lines) appear as straight lines 
in this case. The different slopes are due to the different nonlin- 
ear structure growth (i.e., effects of the second order term on the 
BAO) depending on cosmology. The predictions for T256 (short- 
dashed magenta) and G576 (long-dashed black) overlap. Note that 
a fair fraction of the difference between SSEW and the other sets 
is explained by the difference in cosmologies. 



are almost the same. Again, G576 and T256 follow a 
D^(z) dependence at high redshifts (i.e., at small values 
of D(z)^) while the dependence becomes weaker at low 
redshifts. 

In redshift space, the displacements of mass tracers 
are larger along the li ne of sight direction by (1 -I- /) 
(jEisenstein et al.ll20Q7| ). and we therefore expect more 
shift along the line of sight by the amount that depends 
on /, where / is the logarithmic derivatives of the linear 
growth factor (~ ri„i(z)°-^). From G576, we derive shifts 
in redshift space that are greater than real space by 34%, 
23%, and 17% at z = 3, 1, and 0.3, respectively. This 
roughly agrees with the shift increase of ~ 1+//3 = 33%, 
30%, 25% at z = 3, 1, and 0.3 that would be expected for 
a spherically averaged redshift space power spectrum if 
the shift along the line of sight were (1 + f) times larger. 

4. EVOLUTION OF PROPAGATOR 

The degradation of BAO in the final density fields, 
before reconstruction, can be reasonably approximated 
by a linear correlation function convolved with a Gaus- 
sian function in configuration space or equivalently a 
Gaussia n damping of th e linear BAO peaks in Fourier 
space (lEisenstein et all [200l IC rocce fc Scoccimarrd 



2b06VEiscnstcin et al."'2007';'Crocce & Scocciniarro|l2008t 
Matsubara 2008: Padmanabhan ct al. 2009), based on 
the Zel'dovich approximation. Meanwhile, it has not 
been investigated whether such an approximation is a 
good descri ption for the density fields after reconstruc - 
tion (but see lPadmanabhan et al.ll2009HNoh et al.ll2009( ). 
We can assess the exact amount of the remaining BAO 
signal in the final density fields in a numerical way by de- 
riving the cross-correlation betwe en the initial and final 
density fiel ds, i.e., "propagator" (jCrocce fc Scoccimarrd 
[lMl2008h . In this section, we derive propagators (eq. 
[5]) for real and redshift space and before and after re- 
construction; 



G(fc,z) 



< (5iin(fc, z)(5(fc',z) > 

< (5lin(fc, z)(5iin(fc',z) >' 



(5) 



where 5\ia{k,z) is the initial linear fields that is linearly 

scaled to z, and 5{k, z) is from the final density fields at 
z. This implies that the nonlinear power sp ectrum can 
be modeled as (jCrocce fc Scoccimarrd[2008f) 

P,,,{k,z) = G{k,zfPUk,z) + P^c{k,z), (6) 

where Pun is the input power spectrum that is linearly 
scaled to z, and Pmc is the mode-coupling term that 
describes a portion of power spectrum that is not directly 
correlated with the initial fields. 

Figure [7] shows G{k) in various cases. The black solid 
lines and the dotted lines, respectively, show the propa- 
gators from our simulations and the Gaussian damping 
model based on the analyti c Zel'dovich approximation 
from fEiscnstc in et al.l (|20'07[ f|. In redshift space (right 
panels), unlike in real space (left panels), we use a by-eye 
estimate of an isotropic Eni to represent an anisotropic 
damping model. Before reconstruction, based on the 
real-space G{k), the Zel'dovich approximation is a good 

* We derive the estimates of nonlinear damping parameter S„i 
for the pairs separated by the sound horizon scale. 
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Fig. 7. — G{k) at z=3.0, 1.0, and 0.3 in real (left) and redshift space(right), before (black for G576 and red for G1024) and after 
reconstruction (blue for G576 and magenta for G1024). The Gaussian m odel of G(fc) is over-plo tted with dotted lines: Ejji for real space 
before reconstruction is chosen based on the Zel'dovich approximation in lEisenstein et al.l l|20071 ). 



description of G{k) for the region of 0.1 < G{k) < 1, 
except for z = 3 . Meanwhile, it is difficuh to estimate 
the amount of the remaining Gaussian damping after re- 
construction, as the effectiveness of the reconstruction 
will vary depending on redshifts as well as shot noise: 
for example, at high redshift, there is less nonlinearity 
from which to recover BAO. The blue solid lines in the 
figure show the propagators from our simulations after 
reconstruction, and blue dotted lines show our crude by- 
eye estimation of the corresponding Gaussian damping 
model. Although the Gaussian damping model appears 
to be a worse description for G{k) after reconstruction, 
we estimate that the decrease in the nonlinear damping 
parameter due to reconstruction is roughly a factor of 2 
at z = 0.3. At z = 1, we show both G576 and G1024. 



Note that the propagators are very similar before recon- 
struction despite the different force, mass, and volume 
resolution: it differs only by 3 4% at = 0.2h Mpc~^. 
The difference in the propagators is bigger after recon- 
struction: the reconstruction seems somewhat more ef- 
fective in G1024 in real space, probably due to its slightly 
higher amplitude on quasilinear scale and smaller shot 
noise (Figure [T|). 

5. THE MODE-COUPLING TERM, Pmc 

In the previous section, we investigated the propagator, 
which is a description of the signal in the signal to noise 
ratio associated with the acoustic scale measurement. In 
this section, we investigate the source of the observed 
shift of the acoustic scale shown in Figure [2j Assum- 
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show G576 in redshift space before and after reconstruction. The red line shows the observed nonlinear BAO, Pnl,osc/^'lin,smi and the 
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^nl,osc/^liii,sm " G'gj„Piin,osc/^'liii,sm > therefore the source of the observed nonlinear shift. Note that the oscillatory, non-smooth feature in 
the black lines indeed increases with decreasing redshift. The blue line is for the non-smooth contribution from the mode-coupling term, 
i.e. fMC,osc/^lin,smi and the green is the non-smooth contribution from the propagator, 2Gs„-|(Gosc/Gsm)- The features in the blue lines 
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deviation of components from zero on small wavcnumbors in the reconstructed cases is due to a poor fitting to the smooth components. 
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ing that the propagator is smooth, the source of the ob- 
served shift should reside in Pni{k, z) — G{k, z)^Piin(fc, z), 
and therefore the mode-couphng term PmCj by its def- 
inition in equation ([6|). Here P^iik^z) is the measured 
nonhnear power spectrum with its acoustic scale shifted 
relative to the linear power spectrum, and G{k, z)'^Piia 
is a nonlinear power spectrum with the degraded BAO 
but without the shift. Therefore the difference between 
these two power spectra will show the source of the shift. 
More specifically, we are interested in oscillatory features 
in Pni{k, z) — G{k, z)^Piin(fc, z) and Pmc that are out of 
phase in wavenumber relative to the BAO and there- 
fore mimic the derivative of the BAO with respect to 
wavenumber. 

As we are interested in the oscillatory components only, 
we rewrite equation ([6|) and separate the smooth and 
non-smooth components in Pni, G{k) (if any), and Pmc- 
Then 



nLsin 
2 



Gsm 

^'mCosc 



X Plin.sm I 1 



lin.osc 



Pin, 



P 



MC.sm 



(7) 



(8) 



where Pni,sm, flin.sm, Gsm, and Pivicsm are the smooth 
components and Pni.osc, fiin.osc, Gosc, and Pmc,osc are 
any non-smooth components that include oscillatory 
components as well as random noise. Note that we have 
not assumed that Ggm is smooth. We remove all addi- 
tive terms that involve only smooth components, as these 
can be marginalized over, and then keep only the first or- 
der term in Gosc /Gsm- Finally, we divide all terms with 



-Pnl,o 



lin.sm 



Py, 



lin.osc 
lin.sm 



2Gi 



Go 



MC.osc 



P 



lin.sm 



(9) 



Therefore the nonlinear shift that is imprinted in p"''"""^ 

-Oin.sm 



G 



2 Pn 



will depend on 



PmC,: 



— , as long as Go 



0. 



Our goal is first to show the source of the observed non- 
fni.oac ^2 -Phn.oac ^^^j sBCOud, to comparc 



linear shift in 



-G 



sm p.. 



it with 



after inspecting 2Gs„jp^ for any obvious 



oscillatory feature. 
In detail, we derive the smooth components, Pni.sm, 
Gsm, and PMC.sm in the following ways. For 
-Pni,sm, we use equation ([3]) but with Pm^Pri w wher e Pnw 
is the nowigg le form from lEisenstein fc Hul (|1998D . We 
then recycle the best fit a, B{k), and A{k) derived for Pni 
to produce Pni,sm because directly fitting to Pni,sm, which 
does not have any BAO feature, produces degeneracies 
among the fitting parameters. We follow the same proce- 
dure for Piin,sm using the randomized power spectrum at 
z = 50. In order to generate Gsm, we combine a Gaussian 
damping model with a Fade approximant and polynomi- 
alf|3. Finally we derive PMC.sm = -Pni.sm - G^^Pun.sm- 
The resulting smooth fits may have small residuals on 

where a'-s, b'-s, c's, and d'-s are fitting parameters. 



small and large wavenumbers due to the broadband dif- 
ference between Pnw and the original P™. However, we 
are not correcting for this, as we are only interested in 
a qualitative rather than quantitative understanding of 
the shift and any defect in our models of smooth com- 
ponents will reveal itself as a broadband deviation of 
the non-smooth components from zero. Note that by 
the definition of Pmc, non-smooth components or a de- 
fect in the fit to Gsm will be propagated to Pmc,osc and 

Pn,n,r ^2 Plin,, 



G^ 



All these should be considered as 



caveats when we interpret the resulting, non-smooth fea- 
tures. 

The black lines in Figure [8] show the source of the ob- 



served shift of the acoustic scale, i.e. 



G'i 



in real space. The first three panels show G576 in real 
space before reconstruction and the middle-right panel 
shows G576 in real space after reconstruction. The bot- 
tom two panels show G576 in redshift space before and 
after reconstruction. Note that the oscillatory, non- 
smooth feature in the black lines increases with decreas- 
ing redshift, which is consistent with the qualitative 
trend in Figure [2] The green lines show the non-smooth 
component in G(fc), i.e., 2Ggjjj(Gosc/Gsm), and the blue 
lines show the non-smooth component in Pmc, i-e., 
PMC,osc/-Piin,sm. The nou-smooth components of G{k) 
seem much smaller than that o f Pmc , which is consistent 
with IPadmana bhan fc W hit(i ()2009l ). With G(fc) being 
much smoother than other components, Pmc.osc appears 
in good agreement with features in — G'^^ ^im.oac 

and therefore indeed seems responsible for the observed 
shift. 

We do not see a strong trough in Pmc,osc near k = 
0. 057/t Mpc~^ that corresponds to P-node 1 in Figure 6 
of ICrocce fc Scoccimarrol (j2008D . This is probably be- 
cause the actual contribution from this trough is very 
small as the PMC,sm/flin,sm we derived is less than 1/20 
at this wavenumber. Considering the errors on the 
power spectrum (gray error bars) relative to the oscil- 
latory feature and the relatively large contribution to 
the BAO information f rom k = 0.1 — 0.2/i Mpc^^ (see 
ISeo fc EisensteinI (|2003f ) and also see Figure HI] in § 0, 
we estimate that the off-phase oscillatory feature in Pmc 
near k = 0.15/i Mpc~^ substantially contributes to the 
shift. In fact, the rnode-co upling shift measured by 
ICrocce fc Scoccimarrol (|2008[ ) for this node (i.e., P-node 
4 in their Figure 7) is in good agreement with our mea- 
surements of shifts. Features beyond k ^ 0.2h Mpc^"'^ 
are on phase with the BAO, so we do not expect them 
to contribute to the observed shift. The feature near 
k — 0.1ft. Mpc~^ might also be responsible for the shift, 
but given the deviation of Gosc from zero, it is hard to 
judge whether the feature is real or due to an imper- 
fect fit to the smooth component. The middle-right and 
bottom-right panels show the effect of reconstruction at 
z ~ 0.3 in real space and redshift space, respectively. 
Oscillatory features in the black and the blue lines have 
substantially decreased and we no longer find any evident 
off-phase feature relative to the BAO. The broad-band 
deviation of components from zero on small wavenum- 
bers in the reconstructed cases is due to a poor fitting to 
the smooth components. 

Figure [9] shows the result for G1024 a.t z — 1 before 
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Fig. 9. — Non-smooth components of P^i at z = 1 for G575 and G1024. The left panels show G576 and the right panel shows G1024. 
The different small-scale clustering shown in Figure [T] suggests that the overall mode-coupling term is slightly different between G1024 and 
G576. We however do not find an indication for a difference in the non-smooth component of Pmc between G1024 and G576 induced by a 
difference in force and mass resolution or size of the box used for the simulations. 



reconstruction, in comparison to G576. The different 
small-scale clustering shown in Figure [U after account- 
ing for the propagators of G1024 and G576 from § [U 
suggests that the mode-coupling term is also slightly dif- 
ferent between G1024 and G576. An important question 
to ask is whether or not there is a difference in the non- 
smooth component of Pmc between the two sets, as this 
will be an indirect way to test the shift difference between 
G576 and G1024. From Figure IHl there seems to be a re- 
semblance between G576 and G1024 in the oscillatory 
feature of Pmc over k = 0.1 — 0.2h Mpc~^. Meanwhile 
the feature is much less significant in the case of G1024 
relative to its large sample variance (i.e., the gray error 
bars) which is consistent with our null detection of shift. 
We therefore do not find an indication for a difference in 
the non-smooth component of Pmc induced by a differ- 
ence in force and mass resolution or size of the box used 
for the simulations. 



6. SIGNAL TO NOISE RATIO FROM G{K) AND 

Pmc 



6.1. A toy model 



lm(6) 




(=D(z)/D(0) bj 



Re(6) 



Fig. fO. — A schematic diagram of our toy model to show the 
Fourier coefficient of nonlinear density field, Sf, relative to the 
Fourier coefficient of initial, linear density field, for a given 
realization. The projection of Sj onto <5iin is G(fc)5iin and the 
component uncorrelated to (5iin is Smc- 



As pointed out in iPadmanabhan et al.l ()2009D . the re- 
constructed field is not the linear density field at second 
order. In an attempt to visualize the difference between 
the linear field and the reconstructed field, let us intro- 
duce a toy model in which the density field after recon- 
struction or before reconstruction is related to the initial 
density field in the complex plane of Fourier space as 
following: 

6f{k) = G{k)SUk) + SMc{k). (10) 

That is, we describe the final field as a combination of 
a component along the direction of (5iin, G{k)Siim and a 
component that is statistically uncorrelated to the ini- 
tial field, i5mc (Figure fTO|) . Here G{k) is a real-valued 
propagator and Sn^ is a linear density field scaled with a 
linear growth factor, as before. Note that this toy model 
is not strictly correct: for example, we are assuming that 
(5ii„ and ^mc SiVC statistically independent. Nevertheless, 
the toy model above returns a consistent definition of a 
propagator when multiplied with Sim and averaged over 
ensembles: 



< J/(fc)<Si*„(fc) > = (27r)3G(fc)Piin. 



(11) 



In other words, the average of the projection of (5/ onto 

(5iin is G{k)P\in. That is, G{k) ~ 1 means that the average 
projection returns the linear density field. 

In this picture, Pmc represents the random dispersion 
around the point GSnn- 

< {Sf ~ GSlin)iSf - GSlin)* > = < (^MC^MC (12) 
->P/-G2Piin = PMC. (13) 

Therefore, the nonlinear density field, whether before 
or after reconstruction, can be described by G(fc), repre- 
senting the average projection onto the linear field, and 
Pmc 7 representing a random dispersion around the linear 
field, for a single Fourier component. Obviously, in or- 
der to recover the linear density field, we require at least 
G{k) = 1 and Pmc = 0. In the case of BAO, the recovery 
of BAO signal is well-described by G{k) and therefore 
is related to the average projection on the linear field. 
Meanwhile, it is important to note that the recovery of 
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Fig. 11. — Top loft: G{k")^ Pwi^ and Pmc before (thin and thick dotted lines, respectively) and after reconstruction (thin and thick solid 
lines, respectively) in redshift space at various redshift. After reconstruction, G(fc) increases toward unity and Pmc decreases. Top right: 
the corresponding signal-to-noise ratios per Fourier mode in redshift space derived from G(A;)^Piin and Pmc- Bottom: the total BAO 



signal-to- noise ratio as a function of the maximum wavenumber femax. 
to the variances at various redshifts from Table [5] 

the broadband shape of the hnear power spectrum wiU 
be hard to diagnose with G(\i) alone, in the presence of 
Pmc that also contributes to the broadband shape by 
definition, and therefore needs to take a different strat- 
egy than ones that are effective for the recovery of BAO. 

Wc can relate G(fc) and Pmc to the signal-to-noisc ratio 
of the power spectrum. Since wc look for a linear feature 
in the final power spectrum, the signal, Psignai, is 

< (<5/ -(Jmc)!^/ -(5mc)* > = G2(fc) < 4n^n„ >(14) 

^Psignal(fc)=G2(fc)Piin(fc), (15) 

as expected and Pmc contributes to the noise through 

6.2. Signal to Noise ratio 

When assuming a Gaussian error, the signal-to-noise 
ratio of the standard ruler test, relative to the linear 
density field, now can be written as 



G2P, 



1 



S _ 

N " G2Pii„ + Pmc " 1 + PMc/(G2Pii„) ' 



(16) 



See the text. The data points are ratios of the variance at 2 = 50 



which is valid whether or not we believe the toy model 
in the previous section. It is evident that the optimal 
reconstruction will increase G(fc) to unity and decrease 
the dispersion, Pmc, over the scale of interest. 

In the top left panel of Figure [TTl we show G{k)^P\in 
and Pmc before (dotted lines) and after reconstruction 
(solid lines) in redshift space, derived from G576. One 
sees that the reconstruction not only moves the aver- 
ages of the Fourier modes near the corresponding linear 
positions, but also decreases the dispersion around the 
averages. In the top right panel, we show the increase 
in the signal-to-noise ratio per Fourier mode in redshift 
space due to reconstruction. For the linear power spec- 
trum, the ratio will be unity over all scales, implying that 
the signal here does not only incorporate the BAO but 
all information in Pun- 

We then derive a cumulative quadratic sum of the 
signal-to-noise ratio after weighing with the avail- 
able number of Fourier modes at each wavenumber. 
Also, in order to single out the BAO information 
in the signal-to-noise ratio, we weigh the signal in 
Piin at each wavenumber with the Silk damping ef- 
fect, i.e., exp [— (fcPsiik)^'^]/Piin where Psiik is the Silk 
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damping scale, before the quadratic summation (see 
iSeo fc Eisensteinll2007l for a similar operation): i.e., 

N G^Pun + Pmc' ^ ' 

In the bottom panel of Figure [TTl we show the total BAO 
signal-to- noise ratio as a function of fcmax- The ratio 
is normalized to unity for the linear power spectrum at 
^max = 0.4/1 Mpc~^, for convenience. The dependence 
on k for small wavenumbers is steeper than which is 
expected from counting the number of available Fourier 
modes, and this is because exp [— (fci?giik)^'^]/Piin also 
increases with k before the exponential damping domi- 
nates. Finally we compare these ratios with the variance 
of the shift at z = 50 relative to the variances at various 
rcdshifts (from Table [5]) that are measured in § |3] us- 
ing fcmax = 0.4/1 Mpc~^ (data points). In redshift space, 
since a spherically averaged P(k) from N-body simula- 
tions suffers a larger noise by 7%, 9%, and 10% at z=0.3, 
1.0, and 3.0 than when u sing the full anisotropic sig- 
nal (jTakahashi et al.ll2009al ) , we correct for the measured 
ctq 's by this factor before comparison (also see § [5] for 
more explanation) . After the correction, we find less than 
12% discrepancy in errors on shift between the derived 
signal-to-noise ratio and the measured errors on shifts. 

The role of Pmc implies that the signal-to-noise ra- 
tio will be different for the reconstructed redshift-space 
power spectrum, depending on whether or not we per- 
form FoG correction (i.e., FoG compression) during re- 
construction. As a re minder, the reconstruc tion tends 
to stretch FoG further (jEisenstein et al.ll2007D . as shown 
in Figure [T] A simple way to reduce this effect is, to 
compress the spacial dimension along the line of sight 
by a typical peculiar velocity dispersion of halos rela- 
tive to the Hubble expansion, identify clusters using an 
anisotropic friends-of-friends algorithm, and move all the 
particles to their center of mass of the cluster along the 
line of sight, decompress the spacial dimension along the 
line of sight, before conducting reconstruction. We find 
that, for a moderate level of FoG compression, which re- 
covers the level of FoG effect before reconstruction, G{k) 
improves slightly, while the increase in Pmc is more no- 
ticeable: as a result, the signal-to- noise ratio worsens 
slightly (i.e., - 5% effect at fc - 0.2/iMpc-i). For 
more intensive FoG compression to remove most of the 
FoG effect, the signal-to-noise ratio decreases as much as 
~ 15% at fc '--^ 0.2/1 Mpc~^ despite the improvement in 
G{k). Therefore, we conclude that, in terms of the signal- 
to-noise ratio, reconstruction without FoG correction is 
most effective, although it probably requires more care 
in the process of an anisotropic fitting due to the strong 
FoG effect. On the other hand, including a moderate 
level of FoG correction in the reconstructed power spec- 
trum might improve anisotropic fitting at the expense of 
total signal-to-noise ratio. 

7. IMPROVING RECONSTRUCTION: 2LPT AND 
ITERATION 

In the previous section s, we used the reconstruc- 
tion method presented in lEisenstein et al.l ()2007f ) that 
is based on a linear PT continuity equation. The final 
density field we use in the continuity equation to derive 
displacements is, however, nonlinear on small scales, and 



we therefore smoothed gravity on small scales by using a 
Gaussian filter. Although most of the displacem ents are 
due to large-scale flows (jEisenstein et al.l I2007D , which 
is the reason why our fiducial method works well, there 
must be a small amount of nonlinear contributions to the 
total displacements: smoothing the small-scale gravity to 
mimic a linear density field makes one side of the conti- 
nuity equation linear while, optimally, we want to derive 
displacements on the other side both due to linear and 
nonlinear contributions. In this section, we attempt to 
improve the performance of reconstruction by adopting 
a few more operations into our fiducial reconstruction 
method. We try 2LPT and iteration. 

First, we modify the continuity equation based on the 
second-order Lagrangian perturbation theory. That is, 
we de rive the displaceme nt fields from the following equa- 
tions (|Scoccimarrolll998l and references therein): 

q = -V^^'^-^V^^'\ (18) 

where q is the estimated displacement fields of the parti- 
cles from the initial positions, and therefore the reverse 
of the displacements we will apply to the real particles 
and reference particles at low rcdshifts, 0^^^ is the poten- 
tial in the linear continuity equation derived from linear 
density field. In our case, we assume that the Gaussian 
filtered nonlinear density field closely represents the lin- 
ear density field and derive Then 

= (^obs,filtorod (19) 

and 0'^^) is the second-order correction: 

v'^^'=j:^^^^-i^^^- (20) 

i>j 

The gradient and therefore q are evaluated in the coniov- 
ing (Eulerian) coordinates in our method, and therefore 
to be exact we are deriv ing an Eulerian displacement 
field (jBouchet et al.lll995[) . In redshift space, we agam 
use equation ([T5)) . but this time 0*^^^ and t/j^^^ arc derived 
based on the redshift-space density fields. As a minor de- 
tail, the Gaussian smoothing filter and window function 
correction due to GIG are applied only once for 4>^^^ . 

As the nonlinear shifts are already close to zero with 
our fiducial reconstruction method, we focus our atten- 
tion on the improvement in the signal to noise ratio, that 
is, the reduction of damping in G{k). Figure \n\ shows 
the effect of adopting 2LPT into our reconstruction for 
different smoothing lengths, in comparison to our fidu- 
cial Zel'dovich approximation-based method. As G{k) is 
the ratio of power spectra from the same random seed, a 
good portion of the random fluctuations is cancelled in 
calculating G{k). We therefore use only one box of G576 
simulations. We further reduce the noise by re-binning 
all power spectra using A/c = 0.005/i Mpc^^ before gen- 
erating G[k). 

Overall, the gain of adopting 2LPT is not significant 
in this specific method: while it pushes G{k) closer to 
unity on large and quasilinear scales, it is only ~ 2% ef- 
fect on G{k) when the fiher size R = IQh'^ Mpc. The 
figure also shows that a smaller smoothing filter helps re- 
construction on small scales while less effective on large 
scales, probably due to using more nonlinear regions of 
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Fig. 12. — Effects of adopting 2LPT into reconstruction (solid lines), relative to our nominal Zel'dovich approximation (ZA) based 
reconstruction (dashed lines). The left panel shows the difference due to 2LPT in real space at z = 0.3 when we use a Gaussian filter size 
of i? = 10h~^ Mpc and R = 7h~^ Mpc. The right panel shows the effect in redshift space. We use only one 8h~^ Gpc^ box but rebin G{k) 
in Ak = O.OOS/i Mpc~^ in order to decrease the noise in G{k). The dotted line is for G{k) before reconstruction and light gray straight 
line shows G(fc) = 1. 2LPT seems to correct for the small remaining deviation in G{k) from unity on large scales, which is not likely to 
make a noticeable change in a signal-to-noise ratio. Note that as we decrease R, that is, as we use more nonlinear scale information in the 
density field to reconstruct BAO, G{k) overall becomes closer to unity on small scales at the expense of more deviations on large scales: 
2LPT seems to help fix this deviation on large scales. 
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Fig. 13. — Effects of iterative reconstruction. G(k) at z = 0.3 for various number of iterations. "II" corresponds to our fiducial 
reconstruction. A filter size of i? = 10h~^ Mpc is used. We use only one 8/i~^ Gpc'' box but rebin G(fc) in Afe = 0.005/i Mpc~^ in order 
to decrease the noise in G(k). The light gray line is for G{k) = 1, and the dotted black line is G{k) before reconstruction, the colored solid 
lines are for G{k) for a different number of iterations. In the left panel, G(k) for 12, 13, and 14 is above 1, and it is due to the amplitude 
of reconstructed power spectrum being slightly larger than the linear theory. The right panel shows the differences after the different 
amplitude is renormalized. 



the density fields. Adopting 2LPT seems to fix the prob- 
lem. We will further investigate variations in the 2LPT 
implementation in future work. 

Second, we adopt iterative operations into our fiducial 
reconstruction. At each step in detail, we derive dis- 
placement fields based on the density field of real parti- 
cles (Jobs) at that step and displace the real particles, 
while the new density field is derived from these dis- 
placed particles. At the end of the iterations, we dis- 
place the uniformly distributed reference particles by the 
sum of the displacements from all the iterative steps. Fi- 
nally we subtract the reference density field ((Jrof.N) from 
the density field of the real particles (Jobs.N) and con- 
struct the reconstructed density field ((5rcc)- The process 
is schematically shown below: 



btep i: dobs.o — > Oobs.i — > dobs, 2 • ■ 



Step 2: Jrcf.O '^^il^?" (5rcf,N 



Step 3: (5r, 



So 



bs.N 



qN- 



-V<?>JV-i 



5o 



bs,N 



Figure [13] shows G{k) as a function of the number of 
iterative steps: eg., "II" means our fiducial method with 
= 1 in the diagram above. Note that the iterative 
operations raise the power of the density fields beyond 
the linear growth, as evident from the constant deviation 
of G{k) from unity on large scales (left panel). When we 
renormalize G{k) to match unity on large scales (right 
panel), the figure implies that iterative reconstructions, 
at least this specific choice, helps on large scales slightly, 
then worsens a little bit on the intermediate scales, and 
then helps on small scales by ~ 3%. It may help to 
combine 2LPT and the first two steps of iterations (i.e., 
12), but, considering the small gain in G{k), we conclude 
that the extra effort will not be worthwhile. 

8. TESTING THE FISHER MATRIX ERRORS 
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Fisher matrix estimates of the uncertainty in the acous- 
tic scale have been useful for designing future BAO sur- 
veys, in part due to the simplicity relative to the N-body 
simulations. Obviously, it is important to ensure, by 
calibrating the resulting estimates against the A^-body 
study, that the analytic Fisher matrix calculations in- 
clude reasonable prescriptions for the various nonlinear 
effects and that it therefore agrees with the real Universe. 
However, a precise comparison between the Fisher ma- 
trix estimates and the A'^body results has been limited 
by the sample variance, that is, the available number 
of subsamples. SSEW08 reports a discrepancy of less 
than 25% betwee n the F isher matrix estimates based on 
iSeo fc EisensteinI (I2007D and the A^-body resuhs, . The 
maximum discrepancy is somewhat larger than the ex- 
pected 1 — a fluctuation on errors of 11% based on their 
40 subsamples. Also, the errors from the A^-body re- 
sults tend to be smaller than the Fisher matrix estimates, 
which is opposite to our intuition that the N-body errors 
will be larger as there will be nonlinear contribution to 
the covariance of power spectru m which is not included i n 
the Fisher matrix formalism in iSeo fc EisensteinI (|2007[ ) . 
This result might be due to estimating the scatter from 
only 40 simulations. However it was also suspected that 
there could be a discrepancy of > lh~^ Mpc between 
the Zel'dovich approximation for the nonlinear parame- 
ter Em and the actual damping in G{k) in their A^body 
results. However, from Figurc[71 the consistency between 
the two appears better than ~ 0.5/iMpc~^, except for 
z=3, which is not enough to account for the discrepancy 
we observed in SSEW08. 

In this section, we again compare the Fisher matrix es- 
timates with the A-body resuhs, T256. As in SSEW08, 
we measure the average and the scatters among jackknife 
subsamples but with 5000 subsamples (4999/i^^ Gpc"^ 
per subsamplc) this tii ne. This will all o w us to cali- 
brate the formalism in ISeo fc EisensteinI ()2007f ) with a 

sample variance of l/-\/(2 x 5000) ~ 1%, in principle. 
We find small discrepancies among errors from different 
resampling methods, and they are in general within their 
expected uncertainties. 

We use WMAP3 cosmology in deriving the Fisher ma- 
trix estimates and Em from the Zel'dovich approximation 
for the pairs separated by the sound horizon scale. While 
SSEW08 derive Pghot from a numerical shot noise ex- 
pected from the number density of the tracers, we make 
a more conservative choice of Pshot • We derive the con- 
stant shot noise contribution Pshot from Pni — D^{z)P\in 
at = 0.2ft, Mpc~^, where D{z) is the linear growth 
factor: that is, the derived Pshot accounts for both the 
numerical shot noise and the effects of nonlinear growth 
near k = 0.2h Mpc~^. Note that our choice of Pshot will 
increase the Fisher matrix error estimates relative to the 
choice in SSEW08, therefore taking a more conservative 
stand in resolving the discrepancy reported in SSEW08. 
We derive an effective shot noise rieff = 1/Pshot and list 
"■cffPo.2 in Table [3l 

From Table |3l the discrepancies in error estimates are 
within 10% in real space but at the level of 7%, 15%, 17% 
in redshift space at z = 0, 1, 3, respectively. While we do 
not observe a consisten t trend in real spa ce, th e values 
based on the method in ISeo fc"^i senstein ()2007D tend to 
be smaller than the measured redshift-space errors. How- 



ever the measured redshift-space errors may be slightly 
overestimated for the following reason. In the Fisher ma- 
trix method, we assume a fitting to an anisotropic power 
spectrum, derive errors on Da{z) and H{z), and then 
project these two dimensional errors on the monopole 
mode. Meanwhile, in the A-body results, we fit to 
a spherically averaged redshift-space power spectrum. 
iTakahashi et al.l ()2009a|) pointed out that we should ex- 
pect ~ 1.2 times more variance in the monopole power 
in the latter case at z = 3 (i.e., crf,^ =~ 1.2P^/{Nk/2) 
not P^/{Nk/2), see their eq. [9]). That is, fitting to a 
spherically averaged power spectrum in redshift space, 
as we did with the A-body result, does not perform an 
optimal extraction of a two-dimensional information. We 
therefore should expect that the scatter from the A- 
body results in redshift space is about vL05 = 1.02, 
■\/l-16 = 1.08, and \/L2 = 1.1 times more than the 
Fisher matrix estimates at z = 0, 1, and 3, respectively. 
After accounting for this, we end up with less than 10% 
discrepancies in redshift space as well. This result is con- 
sistent with what we find from the comparison between 
the cumulative BAO signal-to-noise and the measured 
errors on shifts, presented in § |6l which is reasonable as 
the Fisher matrix estimates are based on the propagation 
of errors from power spectrum to the acoustic scale. 

Note that, in real observations, we aim to perform a full 
two-dimensional analysis and therefore it is legitimate to 
use the Fisher matrix estimates that come without this 
extra factor. We could further improve the agreements 
by using a Em measured directly from the A-body re- 
sults, rather than using the analytic Zel'dovich approxi- 
mation, as observed in Figure[7l However, the purpose of 
Fisher matrix calculations is to bypass the high cost A- 
body operations, and the Zel'dovich approximation has 
its virtue for being simple and analytic. We therefore 
conclude that the Fisher matrix formalism in SE07 and 
the A-body results are in a good agreement and the dis- 
crepancy is less than 10% due to an approximation for 
Em and for the covariance matrix in the former. 

9. CONCLUSION 

We summarize the results presented in this paper. 

1. We have measured shifts of the acoustic scale due 
to nonlinear growth and redshift distortions using three 
sets of simulations, G576, G1024, and T256. which differ 
in their force, volume, and mass resolution. The mea- 
sured shifts from the various simulations are in agree- 
ment within 1.5f7 of sample variance. We numerically 
find a{z) - 1 = (0.295 ± 0.075)%[D{z)/D{0)]^-^^^°-^^ 
based on G576. If we fix the power index to be 2, as 
expected from the perturbation theory, we find the best 
fit of a{z) - 1 ^ (0.300 ± 0.015)%[i:»(z)/D(0)]2. 

2. We find a strong correlation with a non-unity slope 
between shifts in real space and redshift space. Mean- 
while the correlation with the shifts at the initial redshift 
is weak. Reconstruction not only removes the mean shift 
and reduces errors on the mean, but also tightens these 
correlations. After reconstruction, we recover a slope of 
near unity for the correlation between the shifts in real 
and redshift space, and restore a strong correlation be- 
tween the shifts at the low and the initial redshifts. We 
believe that, as the reconstruction removes shifts due to 
the second-order, nonlinear process in structure growth. 
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the remaining shifts arc dominated by the initial condi- 
tions. 

3. We find that a propagator is well described by the 
Zel'dovich approximation: for z < 1, we find that the 
discrepancy, is less than O.bh^^ Mpc. At high red- 
shift, Zel'dovich approximation seems to slightly under- 
estimate the amount of nonlinear damping. 

4. We have compared our measurements of shifts from 
fitting with an oscillatory feature in Pmc and find a 

qualitative agreement. 

5. We construct the signal-to- noise ratio of the stan- 
dard ruler test based on the measured propagator and 
mode-coupling term assuming a Gaussian error, and find 
a consistency with the measured errors on shift. We point 
out that the propagator describes the average projec- 
tion of the nonlinear density field onto the linear density 
field while the mode-coupling term describes any random 
dispersion uncorrelated to the information of the linear 
density field. In the case of BAO, the recovery of BAO 
signal is well-described by G(fc), therefore, the average 
projection on the linear field, while the mode-coupling 
term contributes to the noise. On the other hand, the 
recovery of the broadband shape of the linear power spec- 
trum will be hard to characterize with G{k) alone, in the 
presence of Pmc- 

6. We have attempted to improve our reconstruction 
method by implementing 2LPT and iterative operations. 
We find only a few % improvement in G{k). We will fur- 
ther investigate variations in the 2LPT implementation 
in future work. 
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TABLE 1 

W-BODY SETS USED IN THIS PAPER. 



Sample 


Box size 


particles 


Force resolution 


Total Volume 


Code 




(ft-i Gpc) 


{h~^ Mpc) 


{h-'' Gpc^) 




G576 


2 


576^ 


0.1736 


504 


Metchnik & Pinto 


G1024 


1 


1024^ 


0.0488 


44 


Metchnik & Pinto 


T256 


1 


256^ 


3.9062 


5000 


Gadget-2 PM 



Note. — Parameters of the simulations used in this paper. 



TABLE 2 
The evolution of q. 



G576 before Reconstruction 




Qr - 1(%) 


ae-l(%) Aasr{%) 


AQr,z50(%) 


Aas,z5o(%) 


0.3 


0.2283 ± 0.0609 


0.2661 ± 0.0820 0.03774 ± 0.0375 


0.2272 ± 0.0550 


0.2650 ± 0.0781 


1.0 


0.1286 ±0.0425 


0.1585 ±0.0611 0.0299 ± 0.0323 


0.1275 ±0.0341 


0.1574 ± 0.0559 


3.0 


0.0435 ± 0.0293 


0.0582 ± 0.0402 0.0147 ± 0.0213 


0.0424 ±0.0175 


0.0571 ± 0.0318 


50.0 


0.001 ±0.0218 








G576 after Reconstruction 


0.3 


-0.0037 ±0.0298 


-0.0015 ±0.0314 0.0022 ± 0.0128 


-0.0048 ±0.0155 


-0.0026 ±0.0194 


1.0 


-0.0037 ± 0.0269 


-0.0024 ± 0.0305 0.0013 ± 0.0128 


-0.0048 ±0.0113 


-0.0035 ± 0.0182 


3.0 


-0.0021 ± 0.0249 


-0.0058 ± 0.0280 -0.0037 ± 0.0125 


-0.0032 ± 0.0074 


-0.0069 ± 0.0145 


G1024 before Reconstruction 


1.0 


-0.112 ±0.163 


0.002 ±0.233 0.114 ±0.113 


-0.142 ±0.144 


-0.028 ±0.226 


50.0 


0.030 ± 0.088 








G1024 after Reconstruction 


1.0 


-0.093 ± 0.094 


-0.055 ±0.101 0.039 ±0.043 


-0.123 ± 0.044 


-0.085 ± 0.068 



Note. — Fitting range: 0.02/i Mpc^^ < < 0.4/i Mpc-^ We use 504/1"^ Gpc^ for G576 and 44/1"^ Gpc^ for 
G1024. A subscript "r" means a value in real space and "s" means a value in redshift space. 



TABLE 3 

The Fisher matrix estimates in comparison to the N-body data from T256. 





A'^body data 


Fisher matrix 




z 


Sample 




<7,(%) 




n-eff-Po.2 




Real space 


0.0 


T256 


8.27 


0.0217 


8.27 


4.38 


0.0234 




1.0 


T256 


5.26 


0.0133 


5.26 


10.6 


0.0136 




3.0 


T256 


2.78 


0.0090 


2.78 


38.0 


0.0081 




20 


T256 


0.0 


0.0069 


0.53 


100 


0.0072 




z 


Sample 




<7c(%) 




noffPo.2(M = 0) 




Redshift space 


0.0 


T256 


10.0 


0.0274 


8.27/11.76 


4.38 


0.0257 




1.0 


T256 


7.0 


0.0192 


5.26/9.55 


10.6 


0.0167 




3.0 


T256 


4.0 


0.0117 


2.78/5.47 


38.0 


0.0100 



Note. — N-body results are derived by using a total 5000 jackknife samples. 
Values of Em in the fourth column represents the nonlinear smoothing used for the 
template Pm(k) in the analysis of the A'^body data, and Eni and Exy/E^ in the sixth 
column are derived from the Zel'dovich approximations and represents the nonlinear 
degradation of the BAO assumed in the Fisher matrix calculations. 



